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Foreward 


This Final Iteport fulfills the obligation of the Polytechnic Institute of 
New York under NASA Grant GAGW-SI. The period of performance was March 1, 1980 
to February 28, 1981. Dr. Stanley H. Gross was Principal Investigator. The 
program was joint with the Smithsonian Astrophysical Observatory (SAO) tmder a 
subcontract from the Polytechnic Institute of New York. Dr. Mario Gross! led 
the SAO effort. SAO was the source of the Doppler data from the Apollo-Soyuz 
mission. 
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Abstract 


The ASTP Doppler data t«re recalibrated, analysed and related to geophysical 
phenomena and found consistent. Spectra vere confuted for data Intervals cover- 
ing each hemisphere. As many as 14 such intervals were analysed. Wave structure 
is seen' in much of the data. The spectra for all those Intervals are very 
sisiilar in a number of respects. They all decrease with frequency, or with de- 
creasing wavelength. . Power law fits are reasonable and spectral indices are 
found to range from about -2.0 to about -3.5. Both large scale (thousands of 
kilometer) and n^dium scale (hundreds of kilometers) waves are evident. These 
spectra are very similar to spectra of in situ measurements of neutrals and 
ionization measured by Atn»sphere Explorer C. There is just one data sample 
from AE-C from the only orbit for which neasurements are available for the day of 
the ASTP flight. Unfortunately, a lack of receiver lock at this tline produced 
no useful ASTP data for comparison. The nearest available data sample, about 50 
minutes later, has a somewhat different spectral index and wave pools than the 
AE-C data. Its spectru. is made up of a different group of outstanding waves 
as well. ISIS-II data were briefly reviewed. 
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1* The Investigation 

The purpose of the program was to study the Doppler data measured during 
the Apollo-Soyuz Test Program (ASTP) (Weiffenbach et al, 1977), on July 24, 1975 
for information on traveling ionospheric disturbances (TID's) and to evaluate 
the experiment as a measuring scheme. The satellites were at an altitude of 
about 250 km, below the F peak in an orbit inclined approximately 50*. The 
Doppler technique, used between two spacecrafts, enables ionospheric measure- 
ments remote from the perturbing Influence of the spacecrafts and their wakes. 

The differential Doppler measures the rate of change of the column content of 
ionization between the spacecrafts. It had also been planned to seek corrobor- 
ative data from ground measurements to confirm the Doppler measurements and to 
deduce wa/e characteristics. Many of the efforts required to fulfill these 
purposes were completed. Soote were not completed, requiring more time for 
searches and investigations. This report summarizes the accomplishments. 

The tasks were divided between the Polytechnic Institute of New York and 
SAG. SAG was concerned with the proper calibration of the Doppler data, their 
reduction and processing, and the evaluation of the technique as an ionospheric 
measuring tool. The Polytechnic Institute of New York was concerned with the 
management of the effort, the interpretation of other measured data, if avail- 
able, for corroboration, the study of ionospheric structure as detected by the 
measurements, and the use of the technique for detecting TlD*s. 

SAG reduced and processed the data and identified a n»imber of well-known 
ionospheric phenonKna, including the equatorial anomaly, day-night transitions, 
medium-scale TID's and other features, providing confidence in the ability of 
I the technique to measure the ionosphere. The results of their efforts were re- 

I ported in a final report to the Polytechnic Institute of New York, which is 

^ attached to and becoiMs a part of this report. 

f Polytechnic's efforts are detailed in the remainder of this report. Initial 

enqthasls was put on simulating the effects of TID's on ionograms for the purpose 
of analyzing real ionograms. This effort was primarily tutorial. Subsequent and 
major emphasis was on spectral analysis of the Doppler data via the MEM process 
(Ulryck and Bishop, 1975) for wave structure and its interpretation. A study 
was also made of Atmosphere Explorer (AE) data for the same day, July 24, 1975. 
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' The only data were from one orbit of AE-C. Some ISIS-II data were also 

f examined. Ground based lonograms were obtained by S^. For the particular 

f day of the ASTP these lonograms were fotmd to be poor and of little use for the 

purpose. 

Spectral Analysis 

Approximately 9 orbits of data were available from ASTP. A spectral 
analysis of all the data was made by SAO and shown in their attached report. 

Taking all the data together can be misleading from the standpoint of geophysics. 
Ideally, one would want to exar iu: data for various latitude bands. Because data 
are available in 10 second samples only, about 500 data points are available per 
orbit. Since one would like to have at least a few hundred data points in making 
a spectral analysis, the best one can do from the standpoint of geophysics is to 
analyze on a hemispheric basis. Thus, the data were divided into sixteen groups, 
each representing a hemisphere, and each consisting of about 256 data points. 

An example of such a sample of data is shown in Figure 1 which shows Doppler 
shift in Hz vs. a normalized scale representing tin». The normalized quantity 
varies from zero to 2, which corresponds to a time interval of 2560 seconds. The 
start of the data is at 22636 seconds (UT). This is for the northern (summer) 
hemisphere excursion of the spacecrafts in the longitude range 45*Eto -150®W. 
Wavelike oscillations are evident in this snmple of Doppler data. These 
oscillations have half periods from about 60 seconds to 200 seconds, corresponding 
to in-orbit distances of about 450 km to 1500 km for a satellite velocity of 
7.8 km/ sec. Smaller scale structure are also evident. The positive to negative 
Doppler extremes are no greater than about 2.5Hz. 

The power spectrum of these data are shown in Figure 2 plotted as dots at 
discrete frequencies that are Integral multiples of the fundamental frequency, 
1/2560=3.906x10 Hz, based on the duration of the record in Figure 1. The log 
of the spectral power is plotted against the Integral number, or harmonic number, 
i plotted on a logarithmic scale. The highest harmonic number is 128 which 

corresponds to the Nyquist frequency for 10 second samples, 0.05Hz. Thus the 
graph contains 128 dots, harmonic numbers or frequencies. The first dot is on 
t the vertical axis and not visible. The first visible dot Is the second harmonic. 
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A straight line is also drawn In the figure. This is a least squares fit to the 
data and demonstrates that a broad ipectriim is involved decreasing with fre- 
quency with a spectral index^AF in the figure for this case, equal to -2.4371. 

The spectrum is seen to oscillate about the fit with a broad spectral range 
starting at harmonic 14 and continuing to harmonic 42 or 43, corresponding to 
the wavelength scale from 487 km to 1426 km. This range agreeo with the 
periodicities seen in Figure 1. 

A somewhat different sample of data is shown in Figure 3 plotted against the 
normalized time scale as in Figure 1. The time interval is the same, 2560 seconds, 
and the start is at 35916 seconds, UT. This sample is for southern hemispheric 
portion of the orbits in the longitude range 169“E to -27*W. Much less perio- 
dicities are evident in the figure than in Figure 1, much more small scale 
structure is evident. The larger periodicities, however, are typical of those 
seen in Figure 1. 

The spectrum for the data in Figure 3 is shown in Figure 4 together with a 
power law fit. The spectral index of the fit AF“- 2.5494. Some peaks are 
observed in the spectra at the harmonics 7, 14, 17, 20, 27 and 71, corresponding 
to wavelengths of 2853 km, 1426 km, 1175 km, 998 km, 740 km and 281 km, respec- 
tively. In view of the nature of the Doppler data in Figure 3, it is not clear 
whether all these peaks are real or due to the peculiarities of the MEM procedure. 
Peaks at harmonics 14, 27 and 71 may be real, though not precise in frequency, 
due to the discrete nature of the spectra and inaccuracies in the MEM spectral 
estimate. There is a suggestion of nulls in the spectrinn at about the 54th 
and 94th harmonic. The power in the spectrum of Figure 4 is con?>arable to that 
of Figure 2. 

Other examples from both hemispheres are quite similar. Some are shown in 
Figures 5 through 8, Figure 5 and 7 are Doppler data. Each example of the data 
is followed by its spectrum in Figures 6 and 8. All data are plotted vs. a 
normalized time scale corresponding to 2560 seconds. The starting times for the 
data of Figures 5 and 7 are in seconds, UT, respectively, 25286 and 27936. 

Since the data samples are finite, the first 4 or 5 frequencies of their 
spectra are not reliable. Tl»e spectra are better estimates at higher frequencies. 
Therefore, power law fits based on all frequencies, as in Figures 2,4,6 and 8 
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are not valid. It can be aeen In the figures that the fits are not representative 
of the spectral spread at the higher frequencies. The fit should be based on 
spectra starting from a higher frequency or harmonic* such as the tenth. The 
spectra were refitted this way. An example for the data of Figure 3 Is shown 
In Figure 9. This Is the same spectrum as In Figure 4, except that It starts 
from the tenth harmonic. The power law fit Is much more representative of the 
spectrum. Its spectral Index AF differs from that In Figure 4 and has a more 
negative value, -2.9972. The variance of the spectrum about the fit Is given by 
EF and Is 4.8064 In this case. Oscillations about the fit are more clearly seen 
In this figure. Two other examples are shown in Figure 10 and 11. The spectral 
index AF*- 2,3408 and variance EF^4.8183 In Figure 10, whereas the corresponding 
values are -2.0565 and 4.8595 for the spectrum In Figure 11. For ten cases 
studied in this fashion, the spectral Index varied from -2.0565 (Figure 11) to 
■3.41067 with an average .of -2.7148. The variance about the fits varied slightly 
from 4.8034 to 4.8721 with an average of 4.8356. The power levels of the spectra 
are within a factor of 2 of each other In most cases. In one case the power 
level Is within a factor of 4 of -he others. The magnetic activity conditions 
for the time of these measurements were relatively quiet. The magnetic index 
Ap ranged from 5 to 7 and from l"*" to 2 for Kp. The AE Index ranged from 
53 to 475 over the period, but it is not clear whether there Is any connection 
with this index. The uncorrected 10.7 cm solar flux F10.7 %«s 77 for the day. 

It Is evident that both large scale and medium scale spectral components are 
present in all the data analyzed. These may be gravity waves, the large scale 
possibly arising from auroral phenomena and the medium scale from local sources. 
The pattern persists, but no particular components stand out more than others In 
all the samples. The energy levels appear to be about the same, regardless of 
hemisphere. 

Atmosphere Explorer 

Only data from one orbit (8108) of AE-C on July 24, 1975 are available, and 
these only during the time from 45900 seconds (UT) to about 47050 seconds (UT). 

The altitude varied from about 303 km to about 317 km, putting it close to the 
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F peak or above, except in equatorial regions. The longitude range la from 
-151*W to -bl’W. The satellite was entirely In the northern hemisphere during 
this time. ASTP Doppler data during two Intervals that overlap the AE data were 
not useful due to a failure of receiver lock. These «wre from 43906 seconds (UT) 
to 46456 seconds (UT) and from 46566 seconds (UT) to 49116 seconds (UT). A 
comparison can only be made in the next interval for which data is useful, from 
49226-51776 seconds (UT). 

Figure 12 is a plot of electron density measured by the Cylindrical Electro- 
static Probe, CEP, (Brace et al, 1973) on AE-C. Electron density relates to the 
Doppler shift, since the latter Is proportional to the rate of change of column 
density of the electrons between the spacecraft. Wave structure is evident in 
the data. Figure 13 is a graph of nitrogen, helium and argon densities as 
measured by the neutral Atmosphere Temperature ExperlTOnt, NATE, (Spencer et al, 
1973) on AE-C. Nitrogen is the upper curve. Helium is in the middle and Argon 
is at the bottom of the graph. There is some wave structure evident in these 
neutral densities. The measured ion and electron temperatures are shown In 
Figure 14. Tlie electron temperature is the upper curve and is measured by CEP. 

Tlie ion temperature is measured by the Retarding Potential Analyzer (RPA) (Hanson 
et al, 1973). Wave variations are very evident in these a.s?asurements. 

Doppler data from ASTP are shown in Figure 15 plotted vs. the normalized time 
scale for the interval 49226 seconds (UT) to 51776 seconds (UT). These data 
are 2550 seconds in length and are for the northern hemisphere in the longitude 
range -68®W to +97®E. Considerable wave structure is evld--*nt. Periodicities of 
the order of 160 seconds and 126 seconds are evident, con -spondlng to wavelengths 
of about 1230 km and 986 km, respectively. Smaller struct ^ire are also evident. 

The spectrum of these data is shown in Figure 16. Peaks c .-ident in the harnwnlc 
range 16 to 24, correspond to the periodicities seen in Figvire 15. The slope of 
the spectral fit for harmonic 10 and above is -2.6322. Otiier higher frequency 
peaks are evident corresponding to wavelengths in the rang’ about 250 km. 

The spectrum of the AE-C measured electron density is shown in Figure 17. 

A considerable manber of peaks are evident. Two lines are dra%m on the figure. 
This graph is a plot of the log of power vs. log frequency in Hz. Other scales 
shown are the frequency, wavelength in kilometers and harnK’’nlc ninnber. This 
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spectrum is for AE-C data spaced roughly l.S seconds apart for 750 seconds 
starting at 46250 seconds (UT) . The spectrum is estimated using the MEM 
procedure. 

The two straight lines are least square fits for different frequency 
intervals. One has a spectral index of -.975 and is obtained for the frequency 
interval from 20.8 mHz to 130.2 mHz. This interval is approximately 80Z of the 
frequency scale. The second fit has a spectr.'>l index of -1.85. It is obtained 
from the frequency Interval from 20.8 ndlz to 65.1 ndlz, which is half the Interval 
of the first, vniich of the two fits is preferable is not certain. In any event 
this spectrum is flatter than that in Figure 16. This difference may be attribut- 
able to the difference In time and the differences in longitude. Peaks are 
evident in Figure 17 at wavelengths of 599 km, 300 km, 214 km, 146 km, 122 km, 
etc. These wavelengths seem to bear no relationship to those in Figure 16. 

Again the separation in tio» and longitude may be Important. However, it is 
significant that the Doppler spectra and those of AE-C are similar in general 
features. 

ISIS-IT Data 

ISIS-II (G.G. Shepherd et al, 1973) data for July 24, 1975 data were examined. 
Several orbits overlap in time with ASTP from orbit 19950 to 19958. The satellite 
was at an altitude of about 1400 km. lonoprams are of the topside of the ionos- 
phere. Thouth ASTP was below the F-peak, it is of interest to relate the topside. 
Spread-F conditions indicative of structure arc evident in the ionograms. The 
profiles were not reduced, and the record was only briefly examined. No positive 
connections were possible with such limited effort. Itorc analysis is needed. 

Staff 

Computer facilities of the Polytechnic Institute of Nc^ York were used 
during the program. One student was engaged in progransBing. and analysis as part 
of a project effort. 
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Papers 


A paper by R.D. Estes and M.D.Grossi of SAO has been prepared baaed on SAO 
«K>rk. It has been submitted to GFX for publication. Another paper by S.H.GrLSSt 
C. Friedman and M.D.Grossi is in preparation based on tlM Polytechnic spectrdl 
%K>rk. 

* 

Reconmenda t ions 

The program should be continued to seek corroborati^ra data from ground 
facilities to deduce the wave characteristics. Iliis teclmique should be more 
fully utilized for ionospheric measurements when two or more spacecraft are 
available, as may be possible with Shuttle launched equipment. 
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Doppler shift 
22616 seconds 






POWER SPECTRUM. 










The spectrum of the data in Figure 3. All other items 
Figure 2, but for the applicable data. 
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Spectrum of data In Figure 5. All other items are as in Figure 
but for the applicable data 



Another example of Doppler data as In F: 
latitude time Information are contained 





POWER SPECTRUM. 



the applicable data. 
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Figure 9. Spectrum and fit from the tenth harmonic to the Nyqulst frequency 
for the spectrum In Figure A. The straight line fits the data 
better than that In Figure 4. 


POWER SPECTRUM. 



Figure 10. Another spectrum as In Figure 
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Electron temperature (upper graph) and Ion tamparatnra Choc tom eurvm) 
for the same AE-C orbit as in Figure 12. 






Doppler data for a tine pe 
i. The longitude, latitude 








LONG. EXT. 292. 492- >97. 0014 
LAT. EXT. .863->5.047 
FRAM-UORD 104-26/107-62 
U.T. 4.92260E+04/5.17760E+04SECS 
ORDER# = 85 
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